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E E Strong evidence for a gravitational-wave background (GWB) has been reported in the 
nano-Hertz band. Interpreting the origin of this background to be scalar-induced gravi- 
e tational waves (SIGWs), we explore the equation of state (EoS) of the early universe by 
` £ = performing Bayes parameter inferences across the big-bang nucleosynthesis (BBN), cosmic 
microwave background (CMB), and pulsar timing array (PTA) joint observations for the 
first time. Assuming a monochromatic power spectrum for primordial curvature pertur- 
bations, we obtain the spectral amplitude A ~ 1073 — 107! and spectral peak frequency 
f. ~ 1077 — 1075 Hz. We find that the radiation domination with EoS w = 1/3 is com- 
patible with the current observational data, the kination domination with EoS w — 1 is not 
forbidden, while the early matter domination with EoS w — 0 is excluded at more than 2c 


confidence level. These results can be tested with future observations. 
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I. INTRODUCTION 


The equation of states (EoS) is a critical parameter characterizing the thermal evolution process 
of the universe [1]. However, little is known about EoS of the early universe during the era between 
the inflation and radiation domination epochs. Assuming an adiabatic fluid, the density o and 
pressure P can be obtained via the Boltzmann approach [2]. In this context, it can be inferred 
that in the non-relativistic approximation, i.e., p? « E?, the EoS parameter is w — P/p — 0, and 
in the relativistic limit, i.e., E? c p?, we have w = 1/3. Therefore, in practical scenarios, it is 
anticipated that ordinary matter should yield an EoS parameter w ranging from 0 to 1/3. For 
example, in the case of thermal plasma exhibiting strong interaction of QCD, the evolution of the 
relativistic degrees of freedom results in the variation of w within such a range [3-5]. If the early 
universe is filled with a scalar matter field y, the EoS parameter is primarily determined by the 
effective potential V(y) of the scalar field [6-8]. Recently, the oscillating scalar field model has 
been employed to resolve the Hubble tension [9, 10], as well as to address the small-scale crisis 
of large-scale structure [11, 12]. In addition, an epoch between the end of inflation and the onset 
of radiation domination may be early-matter dominated due to primordial black holes (PBHs), 
oscillons, and so on [13-17]. For the early universe, it is reasonable to anticipate that EoS may 
exhibit a mixed behavior, as described in the aforementioned scenarios. It is also well-motivated 


to explore the value of the EoS parameter w through phenomenological studies [18-21]. 


Since the universe is transparent to gravitational waves [22], EoS of the early universe can 
be directly probed via gravitational waves produced in the early universe. It is known that the 
scalar-induced gravitational waves (SIGWs) could be nonlinearly produced by linear cosmological 
curvature perturbations when the latter reentered the Hubble horizon after the end of inflation 
[23-28]. In middle 2023, the pulsar timing array (PTA) data releases showed significant evidence 
for a gravitational-wave background (GWB) in the nano-Hertz band [29-32]. This signal has been 
interpreted as SIGWSs by the authors of Refs. [33—50]. Since the energy-density fraction spectrum 
of SIGWS is sensitive to the parameter of EoS of the early universe [18-20], if such an interpretation 
could be confirmed in the future, we would anticipate that SIGWs are a powerful probe of the early 
universe, particularly for EoS (see review in Ref. [51]). 

In this work, for the first time, we will infer the EoS parameter of the early universe during an 
era between the inflation and radiation domination epochs by conducting Bayesian analysis for the 
North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 15-year dataset 


[31]. Other parameters of the model will also be inferred simultaneously. Comparing our results 


with the existing ones in the literature [52, 53], we will demonstrate the significant impacts of the 
EoS parameter on the inference of these other parameters. On the other hand, once produced 
in the early universe, SIGWs behaved like additional relativistic ingredients, thereby increasing 
the expansion rate of the early universe. Therefore, we take into account the upper bounds on 
the effective number of relativistic species in the early universe from observations of big-bang 
nucleosynthesis (BBN) [54] and cosmic microwave background (CMB) [55]. Analyzing such a joint 
dataset will further refine the above results by more tightly constraining the parameter space. 
The remainder of the paper is arranged as follows. In Section ll, we develop a theory of 
SIGWS considering different EoS parameters of the early universe. In Section IlI, we describe the 
methodology of data analysis used in this work, and present the results of parameter inferences as 


well as their physical implications. In Section IV, we provide concluding remarks. 


II. SIGWS PRODUCED IN THE EARLY UNIVERSE WITH DIFFERENT EOS 


Limited by our current observational capabilities on such small scales in cosmology, we have 
little knowledge of the state of matter in the early universe. In this paper, by assuming an EoS 
P = wp with varying values of w € [0,1] to phenomenologically describe the early-time universe 
preceding the radiation-dominated epoch, we study SIGWSs produced both in the early-time epoch 
and radiation-dominated epoch, as well as the transition from the former to the latter. 

The perturbed spatially-flat Friedmann-Robertson-Walker (FRW) metric in conformal coordi- 


nates is given by 
1 
i= Pe) f- +26) dt? + |1 — 29) &; + shi dadai] (1) 


where ņ is the conformal time, a(7) is the scale factor, ó and w are the Newtonian potential and 
curvature perturbations, respectively, known as the linear scalar metric perturbations, and hij 
describes the second-order tensor metric perturbations, i.e., SIGWs. Due to absence of anisotropic 
stress, we have v ~ $. In the universe filled with perfect fluids with the EoS parameter w, the 


equation of motion of SIGWS is given by 
hy; + 2Hh;, — Ahi = 4A Sap ; (2) 


where Ag? is the transverse-traceless operator, H is the conformal Hubble parameter, and the 
source term Sap consisting of the linear scalar perturbations ¢ is given by 


2(5 + 3w) 4 , a 4 Ko 
304 w) Carer? F 3 (04,0050 Br ah ox) ) T 3 H2 Oud Ong s (3) 


Sab = (1+ w)H (1+ w) 


The evolution of $ is determined by the master equation as follows 
9" E 3(1- w)o' -wAo = 0, (4) 


where we have used the adiabatic speed of sound c2 = w, and focus on 0 < w < 1 for phenomeno- 
logical studies. 

We follow the scenario in Ref. [28] for an extensive study of SIGWs during the transition from 
the early-time epoch of the universe to the radiation domination. The scale factor changes in 
different epochs and can be given by 
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where 5g is the transition time. This also leads to a transition of Hubble parameter from H = 


1 


2[(1 + 3w)n|~! in the early-time epoch to H = 7! in the radiation-dominated epoch. Therefore, 


upon the above setup, we obtain the solutions of h;; in momentum space based on Eq. (2), i.e., 
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where Y,;(x) and J,(r) are Bessel functions of second and first kind, respectively, 8 = —3(1 — 


w)/[2(1 + 3w)] and the superscripts (ET) and (RD) represent “early-time” and “radiation domi- 
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nation”, respectively. Here, we require continuity and differentiability of h;;, ie. h 


and RED = REDY 


ijk T Rijk at NR, such that the time-independent quantities hij can be 


determined as 
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The expressions for the source terms in Eqs. (6) and (7), denoted as se and Se are pro- 


vided in Eq. (3) with ge”) and eR ) defined for the early-time and radiation dominated epochs, 


respectively. Upon solving Eq. (4), we can derive their corresponding expressions, namely, 
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where 9; stands for a random variable related to the primordial curvature perturbations, and oF 


is the confluent hypergeometric function. The continuity and differentiability of øk determine the 


expressions of oe (RD) and ge) as follows 
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It is found that the early-time epoch leaves an imprint on the evolution of linear scalar perturbations 


during the radiation-dominated epoch, which subsequently affects the generation of gravitational 
waves. This suggests that SIGWs in the context of an early-time epoch cannot be simplistically 
understood as simply a superposition of waves across all epochs. The additional influence of scalar 
perturbations on gravitational waves underscores the non-linear characteristics of Einstein's theory 
of gravity. 

The energy-density fraction spectrum of SIGWS at ky > 1 during radiation dominated epoch 
is defined as [23-28] 
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where Pj (k, n) and Pj (k,1)) are defined with the two-point correlations of hij and hi; respectively, 
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In the small-scale limit, where kn > 1, it is permissible to use Ph(k) = Pp(k). The quantities 
/ 

sore and ie ) in Eqs. (15) and (16) are determined by Eq. (9). The stochastic nature of i 

is attributed to that of P. The two-point correlations of 9; can be formally expressed as 
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where P¢(k) is the power spectrum of primordial curvature perturbations ¢. For simplicity, we 
adopt a monochromatic spectrum, i.e., P; = Ak,d(k — ky). Such a spectrum might be related the 
formation of PBHs [56], and has been extensively studied in literature (see review in Refs. [57, 58] 


and references therein). 
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FIG. 1: Energy-density fraction spectra for selected EoS parameters w and transition time rg 


(multiplied with k, to let the combination k,7R be dimensionless). 


In Fig. 1, we depict the energy-density fraction spectrum for various parameters. For the EoS 
parameter w = 1/3, the outcome corresponds to SIGWs in a purely radiation-dominated epoch 
(24, 59]. For the EoS parameter w = 0, the energy-density fraction spectrum experiences an 
enhancement due to the transition to a radiation-dominated epoch. This enhancement mechanism 
pertaining to the early-time matter domination model has also been discussed in Ref. [60]. 

To conduct parameter inferences in the following section, we consider the energy-density fraction 


spectrum of SIGWs in the present universe, which is given by 
h?Qaw.o(k) =~ hOr o x Qaw(k, n) , (18) 


where the expression of Qg@w(k,7) has been given in Eq. (14), and the physical energy-density 
fraction of radiations in the present universe is given as h?Q,9 c 4.2 x 1075 with the dimensionless 


Hubble constant h — 0.6766 measured by Planck satellite [61]. 
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FIG. 2: Posteriors of three independent parameters inferred from the NANOGrav 15-year PTA 


dataset. Here, we consider scenarios with k,rg = 10, 10?, 10°. 


III. PARAMETER INFERENCES FROM THE PTA, BBN AND CMB DATASETS 


In this study, we consider two distinct combinations of datasets. The first dataset constitutes 
solely of NANOGrav 15-year PTA data [31], whereas the second dataset further incorporates the 
BBN [54] and CMB [55] constraints on the effective number of relativistic species. In terms of 
the integrated energy-density fraction, defined as i „dlk h?Naw olk), the upper limits are 
1.3 x 1078 for BBN [54] and 2.9 x 107* for CMB [55]. The lower bound of the integral is denoted 
aS kmin = 2T fmin, where fmin is 1.5 x 10-11Hz for BBN and 3 x 10-1” for CMB [2]. The BBN 
and CMB data have also been taken into account in studies of SIGWs [44] and phase-transition 
gravitational waves [62]. For the energy-density fraction spectrum of SIGWs, as illustrated in 


Eq. (18), we examine three scenarios with the parameter k,7R set to 10, 102, and 10?, respectively. 
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FIG. 3: Posteriors of three independent parameters inferred from the joint BBN, CMB, and 
NANOGrav 15-year PTA dataset. Here, we consider scenarios with k,np = 10, 10?, 10?. 


Consequently, the parameter space under investigation is defined by logi A, logjo( f. /Hz), and 
w. Throughout this work, we take ką = 27 fx. We apply uniform priors for these parameters, 
i.e., logig A € [-3.5,0], log;9(f./Hz) € [-9, —4], and w € [0,1]. Guided by Ref. [52], we perform 
Bayesian parameter inferences across this parameter space, considering both dataset combinations. 
'The BBN and CMB constraints are incorporated by assigning negative infinity to the log-likelihood 


if the integrated energy-density fraction exceeds the specified upper limits. 


By performing Bayes analysis, we show the results for the posteriors of independent parameters 
in Figs. 2 and 3. Here, we have considered the NANOGrav 15-year PTA data for the former, while 
further incorporated the BBN and CMB bounds for the latter. Other settings are the same for the 


scenarios considered in this work. 


When we take into account the PTA data only, based on Fig. 2, we find that the parameter 


A is bounded from lower, but not from upper, i.e., A > 107%. Meanwhile, the peak frequency f. 
is bounded from ~ 1077 Hz to ~ 107? Hz, indicating a micro-Hertz frequency band. Reminding 
that PTA is sensitive to the nano-Hertz band, i.e., ~ 107? — 1077 Hz, and considering the spectral 
profiles in Fig. 1, we find that only the infrared tails of the energy-density fraction spectra can fit 
well the NANOGrav 15-year dataset. Depending on scenarios, loose constraints have been shown 
for the EoS parameter w. 

In contrast, incorporating the BBN and CMB bounds have refined the above results by more 
tightly constraining the parameter space, as shown in Fig. 3. The parameter region, which makes 
the integrated energy-density fraction to exceed the specified upper limits, has been excluded in 
our data analysis. For all scenarios, therefore, we find that the spectral amplitude A is not only 
bounded from lower, but also from upper, i.e., A ~ 107? — 1071. The peak frequency f. is also 
more tightly constrained, i.e., f, ~ 1077 — 1076 Hz. In addition, we still obtain loose constraints 
on the EoS parameter w. However, the most important result may be concerned to exclusions of 
the early matter domination with w = 0 at a more-than-2c confidence level. 

As shown in Figs. 2 and 3, the radiation domination with the EoS parameter w = 1/3 is 
compatible with our dataset combinations for the scenarios considered in this work, and an epoch 


of kination domination with w — 1 [63] is not forbidden. 


IV. CONCLUSIONS AND DISCUSSIONS 


In this work, we conducted a study of the equation of state of the early universe, by assuming 
the SIGW interpretation of the recent PTA data releases. We computed the energy-density fraction 
spectrum of SIGWs during the transition process from an arbitrary w € [0,1] to w = 1/3. See 
Fig. | for illustrating examples. To perform Bayes parameter inferences, besides the NANOGrav 
15-year PTA dataset, we analyzed the dataset combination incorporating this dataset with the 
BBN and CMB bounds. To our knowledge, it is the first time to perform Bayes data analysis for 
the considered topic. T'he results for the posteriors of model parameters have been shown in Figs. 2 
and 3. For the energy-density fraction spectrum of SIGWs, we found the parameter region, i.e., 
the spectral amplitude A ~ 107? — 107! and spectral peak frequency f. ~ 1077 — 10-9 Hz, allowed 
by the joint dataset combination. In particular, we found that the early matter domination with 
the EoS parameter w = 0 is excluded at the > 2e confidence level. However, we showed that the 
radiation domination with w — 1/3 is still compatible with our joint analysis. In addition, we found 


that the kination domination with w — 1 for the early universe is not forbidden by the current 


10 


dataset combinations. The results of this work can be further tested with future observations. 

We generalized the theory of SIGWS to characterize a transition from the early-time epoch with 
an arbitrary w to the radiation-dominated epoch. In fact, the production of SIGWs in the early 
universe with an arbitrary w has been studied in Refs. [18-20]. In these models, gravitational 
waves propagate freely during subsequent radiation-dominated epoch, but are no longer produced. 
However, for SIGWs produced during the early matter-dominated epoch and subsequent radiation- 
dominated epoch, there could be an enhanced production of SIGWs due to a sudden transition from 
the early-matter domination to the radiation domination [60]. It has been shown that the energy- 
density fraction spectrum of SIGWs produced via this enhancement mechanism can dominate 
the total spectrum. Therefore, to study SIGWs in a comprehensive way, we considered SIGWs 
produced in both the early-time epoch with an arbitrary w € [0,1] and subsequent radiation- 
dominated epoch with 1/3, following the approach of Refs. [28, 60]. 

Our findings revealed that the EoS parameter w = 0 is inconsistent with the joint dataset 
examined in this study. This result also showed the importance of BBN and CMB bounds for studies 
of SIGWs produced in the early universe with varying values of EoS. Given that the parameter 
interval of f, is deduced to be roughly 1077 — 1076 Hz, the aforementioned result suggested that 
an extended period of early-matter domination, lasting in excess of 107 seconds (thereby yielding 
k.rg >> 1), is not corroborated by the data. In actuality, when the universe was approximately 
10" seconds old, its temperature would be around 107? MeV. This temperature is lower than 
that which prevailed during BBN [1]. Given the fact that early-matter domination pertains to 
the cosmic evolutionary process spanning from the cessation of inflation to the commencement of 
radiation domination, our research offered fresh perspectives for investigations into early universe 
physics [13-17]. 

The results of this work relied on the assumption that the evidence of GWB in the PTA band 
can be interpreted with the theory of SIGWs. However, such an assumption is not necessarily 
true, since there are also alternative interpretations, as suggested by Refs. [52, 53]. More generally, 
we should consider these possible origins of GWB to make conclusive information of EoS of the 
early universe. In particular, the astrophysical origin of GWB, i.e., the inspiralling supper-massive 
black hole (SMBH) binaries, should be taken into account in the realistic data analysis [53, 64]. 
Considering the above ingredients as well as large uncertainties of the current PTA observations, 
it is challenging to provide confident conclusions to EoS of the early universe at the current stage. 
Nevertheless, our study still showed that SIGWs could be a powerful probe of the early-universe 
physics. 
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